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ABSTRACT 
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identity matri x 


An analysis is presented for multistage multimesh 
gear transmission systems. The analysis predicts the 
overall system dynamics and the transmi ss i bi 1 i ty to the 
gear box or the enclosed structure. The modal synthesis 
approach of the analysis treats the uncoupled lateral/ 
torsional modal characteristics of each stage or compo- 
nent independently. The vibration signature analysis 
evaluates the global dynamics coupling in the system 
The method synthesizes the interaction of each modal' 
component or stage with the nonlinear gear mesh dynamics 
and the modal support geometry characteristics. The 
analysis simulates transient and steady state vibration 
events to determine the resulting torque variations, 
speeds, changes, rotor imbalances and support gear box 
motion excitations. A vibration signature analysis 
scheme examines the overall dynamic characteristics of 
the system, and the individual modal component 
responses. The gear box vibration analysis also exam- 
ines the spectral characteristics of the support system. 

NOMENCLATURE 


CJ] 

<KtiK 
[< s 3 

CKxxl * £C X y] 
[Ky x 3 , [ Kyy ] 

[M] 

R ci 

TF 

Xf.Yf 
x ci *y C i 


rotational (mass moment of inertia) 
matrix 

gear mesh stiffness between i th and k th 
rotor 

shaft stiffness matri x 

bearing direct and cross-coupling stiff- 
ness matrix 

mass-inertia matrix 

radius of gear in the 1 th rotor 

gear generated torque 

gear forces in x- and y-directions 

gear displacements in x- and 
y-directic of the i th rotor 


Ai C t> 


modal function of the i th mode in 
x-direction 


angle of tooth mesh between k* h and i th 
rotor 



modal function of the 1 th mode in 
©-direction 

8j (t) 

modal function of the i th mode in 
y-d 1 rection 

[C xx ],[Cxy] 
CCy x ] , [Cyy ] 

bearing direct and cross-coupling 
damping matrices 

CC t ] 

torsional damping matrix 

F X < t) , Fy< t) 

external excitation forces 

F t (t) 

external excited moment 

F Gx (t).F G y(t) 

gear mesh force in x- and y-directions 

F GT< t) 

gear mesh torque 


[x2],[\fc2] lateral and torsional eigenvalue diago- 
nal matrices 

[$3|(, C<$>t lateral and torsional orthonormal eigen- 
vector matrices of the k th rotor 

INTRODUCTION 

The art and science of analyzing gear transmission 
systems continue to Improve. Power transfer Is neces- 
sary from source to application in mechanical systems. 
Today's engineers and researchers now delve into areas 
of innovative advancement. They seek to quantify, 
establish, and codify methods which can make gear'sys- 
tems meet the ever-widening needs of advancing tech- 
nology. Their objectives nre basic Improvements in 
transmission life, efficiency, maintainability, and 
reliability. They also seek to reduce noise, weight and 
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vibration during transmission operation. Gear transmis- 
sion system studies have Included two main efforts. 

These studies have been on: (1) the localized tooth 

stress/thermal effects during gear interactions, and 
(2) the overall global dynamic behavior of the systems. 

The objective of this paper Is to analyze the over- 
all global dynamics of multistage gear systems using 
localized gear stress/displacement relationships. Equa- 
tions of motion are developed for each gear stage In 
both lateral and torsional directions. Gear mesh force 
and moment relationships couple the lateral and tor- 
sional vibrations and the dynamics of each gear stage. 
Orthonormal modes of the system transform the equations 
of motion to modal coordinates. A self-adaptive varia- 
ble time stepping integration scheme calculates the 
transient dynamics of the system (Choy, 1988). A typi- 
cal three-stage multimesh gear system serves as an 
example. The results of the global dynamics of the sys- 
tem are examined in both time and frequency domains 
using a FFT (Fast Fourier Transform) procedure (Choy, 
1987 and 1988). 

DEVELOPMENT OF EQUATIONS OF MOTION 

For a single stage multimass gear-rotor-beari ng 
system, the equations of motion in the X-Z plane in 
matrix form (Choy, 1987; David, 1987 and 1988) are: 


X F ' Z K tKlKi e c1 - R cK®cK + <x c1 * x cK )cos a KI 
1-1 
\*K 


+ <*d - >cK )s1n a J C0S a Ki (4) 
Summing the forces in the Y-direction results in: 


X F “ Z K tKi[- R ci®li ' R cK®cK + (x ci - X cK )C0S a Ki 
1-1 
UK 

+ l »C\ - *cK )sin a Ki] s1n a Ki <5> 
Summing moment in the Z-direction results In: 
n 

T F * Z R C<[ K tKi ( " R ci e cf - R cK e cK > ] (6) 
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CM3 {x} + CC xx ]{x} ♦ CC xy 3{y} ♦ [K xx + K s ]{x} 


where n is the number of stages in the system. 


+ CK X y ] {y} = {F x ( t) } + {Fq x < t) } 


(1) MODAL ANALYSIS 


In the Y-Z plane, the equations are: 

CM]{y} + CC yx ] {x} + CC yy ]{y} + [K yy + K s ]{y} 

+ CKyxHx} = {F y (t>} + {F Gy < t)} (2) 


To reduce the computational effort, the number of 
degrees-of-freedom of the system is reduced througn 
modal transformat ion . Orthonormal modes for each stage 
result from solving the system homogeneous characteris- 
tic equations. For lateral modes, the equations are: 


Here F x and F y are force excitations from the 

effects of imbalance, shaft bow and support base motion. C M ] { x } + 

^Gx and ^Gy are forces Induced from the gear mesh 
interaction with other stages. 

The torsional equation of motion of the single For torsional modes, the equations are 

st-.ge system is: 


K + K 
xx yy 


{x} 


(7) 


[J][S} ♦ CC T ]{6} ♦ CK T ] {©} = { F T ( t ) } ♦ {F Gt <t)} (3) 

In Eq. 3, F^ represents the externally induced torque 
and F G t represents the gear mesh induced moment. Note 
that Eqs. 1 to 3 repeat for each gear stage. The gear 
mesh forces couple the force equations for each stage to 
each other and the torsional equations to the lateral 
equations (David, 1987; Lin, 1988; Mitchell, 1983). 
Torsional, lateral and interstage coupling relation- 
ships appear in the next section. 

COUPLING IN GEAR MESHES 

Gear mesh forces and moments are functions of the 
relative motion and rotation between the two meshing 
gears and the correspondl ng gear mesh stiffnesses. 

These stiffnesses vary in a repeating nonlinear pattern 
(Cornell, 1981; Pike, 1987; Savage, 1986). The pattern 
repeats with every tooth pair engagement and acts as a 
source of excitation at each mesh. Figure 1 shows the 
coordinate system for the following force and moment 
equations. Summing the forces, which act on the K th 
stage of the system, in the X-directlon results In: 


[J ]{$} + CK T ](9} = 0 (8) 

Averaged X- and Y-direction support stiffnesses bring 
the calculated mode shapes closer to reality. Using the 
modal expansion theorem yields (Choy, 1987): 

{*} - ]T A,^,} (9) 

m 

(y) * 23 MV (10) 

i-i 

( 0 ) * 23 (in 

Inl 

where m is the number of modes used. With the follow- 
ing orthogonality conditions, 
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[♦.]>,][♦.] ■ M 

= [<(> t ] T C J ] * [I] (13) 

the equations of motion in modal coordinates become: 

{A} ♦ Cf 3 T [c xx ] [ 4 >]{A> ♦ C 4 . 3 T [c xy ] C 4 >HB} ♦ [x 2 ] {A} 


m 

' / ^ R cl([ K tKt(' R C) 0 c1 

/ J / U1 

j*1 ( i *K 


' R cK®cK. 


)] 


( 19 ) 


+ [ * )T [ K xy]^HB} - [*] T {F x (t) ♦ F Gx (t» 
{B} ♦ [<t»] T [c yy ][<t,]{B} ♦ Ccf] T [c yx ] E 4 »J{A} ♦ [x 2 ] { 8 } 
♦ [♦^[KyJUHA} = [*] T {F y (t) ♦ F G (t» 

<V ♦ K] T [c T ][« t ]{A t } ♦ [x 2 t ](A t } 

r i T 

- [♦tj {F t <t> + F Gt (t)} 


(14) 


( 15 ) 


(16) 


Thus* the gear mesh force and moment coupling equations 
for the stage in the modal form are: 

m 

E 


[ 4> ] K {X f } - / 

3 = 1 



- ufkj 

[" R ci 6 ci * R cK 6 cK 

L yxj 


where 2 is the station location of the gear in the K th 
stage (Boyd, 1987; August, 1986). 

SOLUTION PROCEDURE 

Rearrange the modal equations of notion developed in 
Eqs . (14) to (16) into: 

X-equat ion 

{A} = -[^[cJuHA} - [*] T [c xy ]c*]{B} - [x 2 ] {A} 

- [*] T [K xy ]t*]{B} ♦ C<J>] T {F j( (t) + F Qx (t)} (20) 
Y-equation 

{ 3 } = -Cct>] T [c yy ] C 4 ,]{S} - [*] T [c y J[*]{A} - [x 2 ] { 8 } 

:A} + [<t>] T {F(t) + F-.(t)} (21) 


1 = 1 
i *K 


and* 0-equation; 

<V— [♦t] T [ c r][*tK>- K]<v 


("cl - *«)“' V * ('cl - » rl ] co - “| 


XI 
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[<t> t ] T {F t (t) + F Gt (t)} (22) 
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U1 

i*K 


+ x 


( x cl - x ck) cos “Kl + ( y ct - y c<) sin “Kf] s1n a i 


XI 


( 18 ) 


A variable time stepping Newmark-Beta integration scheme 
evaluates the modal velocity and displacement at each 
time Interval. In turn* Eqs. 9 through 11 transform the 
modal displacements into absolute/relative displacements 
in fixed coordinates. 

DISCUSSION OF RESULTS 

To demonstrate the application of the analytical 
method, a three-stage multimesh gear system serves as 
an example. Figure 2 shows the geometry of the gear 
system. Stage 1 Is the driver stage. The stage 1 gear 
drives both the output gears directly at a speed of 
1500 rpm. Its input torque of 2.25 kN-m is split 
equally between the two output stages. All the meshes 
are identical with 36 teeth and a contact ratio of 1.6 
as shown In Fig. 3. Although similar, the lateral sup- 
port stiffness for stage 2 is greater than that for 
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stage 3. The system has only minimal Imbalances to make 
the torsional vibration of the system more pronounced 
than the lateral. Figure 4 shows the lateral vibration 
orbits for all three stages at the mesh locations. Note 
that the difference in orbit sizes result from the diff- 
erences In imbalance and shaft stiffness for each stage. 
While stages 2 and 3 maintain a relative circular orbit, 
stage 1 exhibits an elliptical motion due to the tooth 
mesh stiffness interactions with stages 2 and 3. 

Figures 5 and 6 show the gear mesh forces in botn 
time and frequency domains between stages 1-2 and 3-4. 

A substantially larger force is present in the 1-3 mesh 
than in the 1-2 mesh. This is mainly due to the fact 
that stage 3 possesses a higher vibrational magnitude 
than stage 2. The orbi sizes shown in Fig. 4 illus- 
trate this difference. Stage 3's support stiffness is 
less than that of stage 2. 

Table 1 lists both the undamped lateral and tor- 
sional natural frequencies of all three stages with 
their correspondi ng bearing supports. Figure 6 shows 
that the dominating components at stage 1-2 are the 
tooth mesh frequency of 900 Hz and the first stage tor- 
sional natural frequency of 355 Hz. A large average 
force component at zero frequency is also present/ But 
this component does not contribute to the transmission 
vibrations. A considerable force also occurs at the 
input shaft rotational frequency of 25 Hz while a 
smaller component is present in the stage 2 torsional 
mode at 550 Hz. A small force echo is also present in 
stage 2 at the stage 3 torsional natural frequency of 
280 Hz. Similar observations are possible with the 
stage 1-3 force frequency diagram. 

Figure' 7 and 8 show the first three modal contri- 
butions of stage 1 in both the time and frequency 
domains. The correspondi ng mode shapes are In Fig. 9, 
Figure 7 shows both fixed (d.c.) and varying (a.c.) 
responses in each mode. The response magnitude 
decreases with increasing mode number. The frequency 
domain plot of Fig. 8 shows all the modes excited at 
their own natural frequencies. The ratio of fixed 
(d.c.) to varying (a.c.) signal magnitudes in Fiq. 8 
decreases as the mode number increase. Figures i0, 11 
and 12 are plots of the time response, frequency 
response and mode shapes for the first lateral mode of 
each stage. These figures provide the modal information 
of the lateral vibration of all three stages. Note that 
the rotational frequency of 25 Hz is excited in all 
three stages. Figure 11 shows that stage 3, the first 
lateral frequency of 115 Hz is also excited. This is 
due to the characterl Stic of its mode shape which is 
shown in Fig. 12. 

CONCLUSIONS 

This paper presents a vibration signature analysis 
for multistage gear transmissions. The analysis com- 
bines gear mesh dynamics and structural modal analysis 
to study the transmission vibrations. This is a compre- 
hensive method of analyzing multistage gear system with 
effects of geometry complexity, support flexibility, 
mass imbalance and shaft bow. in this method: 

1. The modal method transforms the equations of 
motion into modal coordinates to reduce the degrees-of- 
freedom of the system; 

2. Gear force observations in both the time and 
frequency domains provide good insights Into the source 
of dominating response forces; 


3. Knowledge of modal excitations provide an under- 
standing of the vibrational characteristics of the 
system which can result In Improved transmission per- 
formance and durability; and 

4. The coupling effects of connected structures 
such as the gearbox are easily included in the existing 
modal analysis without sacrificing the above advantages. 
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TABLE I. - SYSTEM NATURAL 
FREQUENCIES 




Mode Stage Stage Stage 

1 2 3 


Torsion natural frequencies, 




550 

280 

610 

820 


Lateral natural frequencies, 



i 


15 

160 

110 

45 

189 

200 

90 

264 

260 



BEARING 
SUPPORTS ^ > 


kth STAGE 


FIGURn 1. - COORDINATE SYSTEM FOR GEAR MESH FORCE AND 
MOMENT. 


i H STAGE III 

FIGURE 2. - TYPICAL THREE STAGE ROTOR -BEARING- GEAR SYSTEM 
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